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[ij Shadowed locations near the lunar poles are almost certainly electrically complex 
regions. At these locations near the terminator, the local solar wind flow's nearly tangential 
to the surface and interacts with large-scale topographic features such as mountains and 
deep large craters, In this work, w'e study the solar wind orographic effects from 
topographic obstructions along a rough lunar surface. On the leeward side of large 
obstructions, plasma voids are formed in the solar wind because of the absorption of 
plasma on the upstream surface of these obstacles. Solar wind plasma expands into such 
voids, producing an ambipolar potential that diverts ion flow into the void region. A 
surface potential is established on these leeward surfaces in order to balance the currents 
from the expansion-limited electron and ion populations. We find that there are regions 
near the leeward wall of the craters and leeward mountain faces where solar wind ions 
cannot access the surface, leaving an electron-rich plasma previously identified as an 
“electron cloud.” In this case, some new current is required to complete the closure for 
current balance at the surface, and vve propose herein that lofted negatively charged dust 
is one possible (nonunique) compensating current source. Given models for both 
ambipolar and surface plasma processes, we consider the electrical environment around 
the large topographic features of the south pole (including Shoemaker crater and the 
highly varied terrain near Nobile crater), as derived from Goldstone radar data. We also 
apply our model to moving and stationary' objects of differing compositions located on the 
surface and consider the impact of the deflected ion flow on possible hydrogen resources 
within the craters. 
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1. Introduction 

[ 2 ] The permanently shadowed lunar craters have a very 
complex environment. These locations are without direct 
sunlight and thus are very cold. The presence of low- 
temperature cold traps located just adjacent to sunlit regions 
allows for the collection of migrating volatiles as they pass 
into darkness, including the distinct possibility' of trapping of 
water and other hydrogen-based compounds [see Vondrak 
and Crider, 2003, and references therein]. We suggest here 
that the polar environment is not only thermally anti chemi- 
cally complex, but also electrically complex with the devel- 
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opment of a plasma miniwake, large surface potentials, 
inward-deflected ion flow and possible dust transport within 
the permanently shadowed craters (PSCs). As we describe, 
the electrical complexity may have a direct effect on the cold 
trap atom collection process. 

[3] The fundamental process described herein is the 
expansion of solar wind plasma into any polar crater. The 
situation is illustrated in Figure !., At the poles, the solar 
wind flows nearly horizontal to the surface and is signif- 
icantly influenced by loca! topography. Specifically, for any 
large mountain or deep crater, solar wind flow will be altered 
(a solar wind orographic effect) whereby the solar wind plasma 
is absorbed on th e upstream side of the obstacle and expands 
into the void formed downstream of the obstructing object 
[ Farrell et ai, 2007J. In essence, a miniwake (i.e., miniature 
plasma wake as opposed to the global-scale wake that forms 
many lunar radii downstream of the moon) is formed on the 
leeward side of such obstructions. 

[4] Plasma wakes/expansion regions are themselves 
complicated [ Crow et a!., 1975; Samir et al . , 1 983; Singh et 
at., 1989; Halekas et al, 2005], Immediately behind the 
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Figure 1. An illustration of the solar wind orographic effect over lunar polar craters. While the solar 
wind flow is initially horizontal, ambipoiar electric fields draw ions into shadowed/solar wind obstructed 
cratered regions. Thus, ion trajectories are deflected into regions that would be considered '‘ion-free.” 
Note the presence of an electron-rich region just leeward of the obstructed flow (called the electron cloud 
region) where electrons migrate into the void faster than the ions. In a steady state situation, the ions may 
never fully catch up to the electrons in these regions. 


obstructing object, a plasma void is established. Thermal 
electrons at the edge of the void (along the flank of this 
void) expand into the region ahead of more massive ions 
because of their higher thermal velocity. Because of this 
charge separation, an ambipoiar electric field (see Figure 1 ) 
develops across the wake flank to retard the electron flow 
and accelerate ions into the crater-generated void. As 
illustrated in Figure i, originally horizontally flowing ions 
will be driven into the crater/void via the diverting force of 
the ambipoiar E field. The ions tend to form beams that have 
a progressive increase in energy to many times the ion sound 
speed, C s [Crow et al, 1975]. As an analogous example, 
such ion beams are commonly detected near the wake flank 
of the space shuttle [Singh et al., 1989] and were in fact 
detected by the Wind spacecraft in the large wake/void 
created antisunward of tine Moon, at a distance of ~7 R L 
[Ogilvie et al ., 1996], 

[<] There are two views of the ion response in the void. 
Samir et al. [ 1 983] presented a self-similar model where the 
expanding plasma remains neutral, i.e., where ions catch up 
with electrons to maintain quasi-neutrality (n, = n,) every- 
where downstream of the obstruction. There is thus charge 
balance in the ambipoiar region. However, if the quasi- 
neutral assumption is relaxed, it is found that directly behind 
the object, the ions never fully catch up, leaving an “electron 
cloud” in the wake region [Crow et al., ! 975], as illustrated 
in Figure I. Particle-in-cell (PIC) plasma codes of the 
plasma expansion process also show this electron-rich 
region [Farrell et al, 2008a] with the “cloud” really con- 
taining energetic electrons from the tail of the solar wind 
electron energy distribution. Compared to the solar wind 
plasma, the electrons in the cloud are velocity -filtered solar 
wind electrons of lower density and in the form of a beam 
at 2-3 times the solar wind electron thermal speed, v to . As 
indicated in the analytical model of Crow et al. [1975] and 
PIC plasma codes {Farrell et al., 2008a], the ambipoiar 
E field reaches a maximum value at the inward limiting edge 
of the ion expansion (called the “ion front”) [see Crow et al, 
1975, Figure 2], The ambipoiar potential is smoothly varying 
throughout the expansion region. 

[6] While the difference between the charge neutral [Samir 
et al, 1983] and nonneutral [Crow et al. , 1975] ambipoiar 


models may appear to be an esoteric academic distinction, in 
reality the model perspective is critical in determining the 
surface potential at the bottom of lunar polar craters and 
behind polar mountains. Specifically, the local surface poten- 
tial is defined by th e current balance of al 1 species reaching the 
surface. In shad o wed/wake regions, there are typically three 
components that primarily define current balance: That from 
the solar wind electrons, solar wind ions, and secondary 
electrons emitted from the surface in response to the incoming 
plasma electrons. In this ease, the shadowed surlace potential, 
ip s , assuming a Maxwellian distribution for electrons and 
ions, a thin sheath, and planar geometry is [Manka, 1973]: 

e<A> ~ kT, lnf j„ ( ! - 6 e ir)/fi) (1) 

where T e is the electron temperature, j Ci , are the electron and 
ion current density, respectively, and is the effective co- 
efficient for secondary electron emission from the lunar sur- 
face that has a value near 0.2-0. 8 for an. electron temperature 
of-lOeV [Halekas et al , 2002, 2009a], The electron and ion 
currents are those located at a point just above the Debye 
sheath along the crater walls and floor. 

[?] Herein lies the dilemma: If the Samir et al. [1983] 
charge neutral model accurately represents ambipoiar pro- 
cesses, electron and ion densities are forced to be equal 
everywhere downstream of the obstacle. The densities be- 
come iie ~ Dj ~ n eo exp(eyLvfcTe) where <p A is the ambipoiar 
potential drop (which is a negative potential). By forcing 
electron and ion concentrations to remain equal, one can 
always obtain current balance between inflowing solar wind 
electrons, surface emitted secondary electrons, and solar 
wind ions, allowing the surface potential to be easily 
quantified via equation (I). In essence, the ions will veer 
from their original horizontal flow and propagate into the 
crater in cadence with the thermally expanding electrons. On 
the other hand, if the Crow et al. [1975] nonneutral ambi- 
poiar model more accurately characterizes wake expansion, 
there will be regions just behind the obstacle that are im- 
mersed in the electron cloud, with no solar wind ions present 
(see Figure 1). An electron-rich plasma will then exist at 
locations where the well-directed ion flow cannot reach the 
surface, especially at locations in the crater where the local 
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Figure 2. The strong negative surface charging within the electron cloud region may be remediated by 
lofted negatively charged dust that acts as a possible current of last resort to keep the potentials from 
becoming extremely large. 


surface normal points away from the flow (i.e., the face of 
the surface is not incident in any way with the flow). As is 
clear from equation (1), absent other current sources, the 
surface potential at such locations goes to negative infinity 
because current balance closure is never achieved. Without 
the ions or some other positive current to the surface to 
balance the electron current, the surface charging goes as 
dQ/dt = J e (1 - 6 e n) exp (eyykT,.) < 0 for all times; the 
surface perpetually charges negative. 

[s] Obviously, this is unrealistic. If the potentials become 
too large, some other near-surface process creates a 
“breakdown” current to provide the required current bal- 
ance. The quantity J; in equation (1) is thus replaced by this 
other current. One possibility is that the secondary electron 
production increases as potentials get progressively larger 
(<5 c f[- approaches unity under larger potentials). This sec- 
ondary production increase may be especially relevant to 
kappa electron distributions, where Te increases progres- 
sively with distance from the wake flank [Halekas et al, 
2005] thereby creating greater secondary production 
(which varies as T e ). In the most extreme case, this sec- 
ondary emission might even take the form of “field emis- 
sion” from surface irregularities. It may be argued that a 
large surface potential at the crater floor will draw in ions 
that are originally flowing at larger distances from the crater 
floor. However, for a crater surface located 10 Debye 
lengths from the ion flow (—150 tn), the surface potential 
will have to overcome a Debye shielding attenuation of 
nearly a factor of 10 6 and for 20 Debye lengths distance, a 
shielding attenuation of nearly a factor of 10 . Only in loca- 
tions close to the ion flow (within a few Debye lengths of the 
ion flow) will there be significant ion diversion to remediate 
the perpetually negative surface charging. 

[»] However, another possible process that has received 
much discussion is lifted/lofted charged dust from the sur- 
face that may act as a natural current source to remediate 
large potential development. Specifically, assuming an 
electron rich plasma reaches the leeward wall of the crater, 
the surface and loose dust will be charged with similar 
polarities (both negative). If surface potentials become large 
and negative, then the local E fields may become large 
enough to loft the like-charged dust [ Stubbs et al , 2006]. This 
lofted dust current J 0 (consisting of upward-moving negative 
charge) may act “as a current of last resort” in place of the 
inflowing positive ions, J ir in the electron cloud region. In 


this work, we will attempt to quantify such a dust current 
and incorporate it into surface charging models. As we 
show, the dust currents should be emitted from the leeward 
surfaces of polar obstacles (mountains, crater walls) where 
inflowing ions are noi directly incident on the surface, 
creating an antisun ward dust flow that may be reminiscent 
of the Apollo 17 Lunar Ejecta and Meteorite (TEAM) 
observations in the Taurus-Littrow valley [Berg et at., 1 976; 
Farrell et al, 2007], The situation is illustrated in Figure 2. 

[to] In this work, we will use topographic relief of the 
lunar south polar region provided by Goldstone radar ob- 
servations [Margot et al., 1999] to quantify the obstacle 
dimensions that creates the solar wind orographic effect. We 
assume that the solar wind plasma has a near-horizontal 
flow' overtop the region, and we then derive values of the 
two potentials (ambipolar and surface potential) that affect 
both electron and ton migration into permanently shadowed 
craters (PSCs) like Shoemaker and Shackleton. We will 
ultimately display a potential map and indicate regions 
where dust lofting and object dust cling may be the greatest. 

2. Topography of the South Polar Region 

[it] Figure 3a shows the Goldstone 3,5 cm (8.6 GHz) 
solar system radar measurements from the south polar re- 
gion of the Moon [Margot et al, 1 999]. The delay in return 
echoes allow a measurement of relative altitude and thus 
provides a topographic map of the region with 600 m res- 
olution. The data are available from the Planetary Data 
System and measurement details are further described by 
Margot et al. [1999]. 

[ 12 ] As indicated in Figure 3a, there is substantial relief in 
the region. The varied topography provides an optimal case 
for developing a solar wind orographic/piasma expansion 
model. To the northeast, near 85S.45E, there are mountains 
that extend to neariy 6 km in height but the terrain also 
appears to be marked by impacts like Nobile crater that 
create abrupt changes in the local relief. Shackleton crater is 
dearly evident in the radar measurement near 89.7°S, 1QQ G 
E. Unfortunately, the radar signal did not map the entire 
floor of this crater. Hence, for crater modeling we use the 
larger Shoemaker crater located at about 88“S, 45°E which 
has a set of clear and distinct echoes from the southern 
portion of the floor. Figure 3b shows a topographic profile 
across the southern portion of Shoemaker crater, with the 
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Figure 3. A topographic map of the south polar region obtained from the Goldstone 8.6 GHz radar sys- 
tem. A profile of the topography across Shoemaker crater is also displayed, and this profile is our primary 
test case used in the study. Regions used to obtain surface potential maps in Figures 7, 8, and 1 0 are 
indicated by blue rectangles. 


crater floor at 4800 ro depth. This particular profile is used 
as our standard test case. We also indicate the locations of 
our regional analysis presented in Figures 7, 8, and 10. 

[ 13 ] We later apply the plasma expansion/surface poten- 
tial algorithm to larger regions, which will prove to be very 
revealing. In order to run such regional models, the deter- 
mination of the ambipolar potential and surface potential 
behind polar mountains and within polar craters has to be 
streamlined. We apply the analytical model of Crow et «/. 
[1975] to map the ambipolar potentiai and we perform a 
very simple surface potential calculation to achieve the results 
at large scale. The process is described further below. 

3. Ambipolar Potential 

[if] Any charged particle reaching the surface down- 
stream from an obstruction has to pass through two negative 
potentials: one associated with the expansion process itself 
(the ambipolar potential) and a second near-surface sheath 
potential defined by local current balance at the surface. Any 
electron has to have enough energy to pass through both 


potentials, and as such, only the most energetic electrons at 
the tail of the solar wind electron energy distribution possess 
the energy required to extend deep behind the leeward sides 
of obstructions. 

[is] As illustrated in Figure 1, the ambipolar potential 
extends from the flank, 'edge into the central region of the 
void. Halekas et al. [2005] measured the large-scale elec- 
tron density' and ambipolar potential across the global lunar 
wake at low altitudes (20-1 15 km) using >6000 passages by 
the Lunar Prospector (LP) spacecraft. They clearly mapped 
out the reduced density and potential structure behind the 
Moon and found that the expansion could be fit to a mod- 
ified self-similar model like that of Samir ei al. [1983], The 
model of Crow et al. [1975] and Samir et al. [1983] have 
similar electron density and ambipolar potential profiles as 
the plasma expands into the void. Without ion measure- 
ments, LP could not determine uniquely if the expansion is 
charge balanced or unbalanced. However, simulations of the 
lunar expansion [ Farrell et al., 1998; Birch and Chapman , 
2001a, 2001b] suggest that there is an imbalance. Both 
cases with be considered herein. 


4 of 14 



E03004 


FARRELL ET AL.: MINIWAKES IN LUNAR CRATERS 


E03004 



Figure 4. The strenglh of the ambipolar potential within 
Shoemaker crater for a westward horizontally flowing solar 
wind. Note that the largest ambipolar potentials exist just 
leeward of the obstacle (i.e., crater’s west wall). The initial 
solar wind conditions that define the properties of the wake 
expansion are n eo = 5/em'\ 1 0 eV temperature, and Vsw ~ 
400 km/s, implying a 15 m Debye length and w pi = 3000/s 
and a 0° elevation angle with the solar wind coming out of 
the west. Table 1 provides more plasma flow details at 
specific points along the crater floor for this particular 
simulation run. 

[is] Both Crow et al. [1975] and recent simulations 
[Farrell et al., 2008a] map out the ambipolar potential (E 
field) behind an obstacle as a function of distance behind the 
obstacle and depth into the void. Both results are compa- 
rable. Crow et al. [1975] provides the calculations in nor- 
malized units that can be easily converted to a general 
relationship assuming a time stationary structure. Fitting to 
those values, we find that the ambipolar potential varies 
approximately as 

ec,\ /k T r = - 1 .08 (V,t ) - I (2) 

where the plasma is expanding into distance t defined by 
Crow et al. [1975] as the normalized depth into a crater, 
Az/Ai> and i is the normalized time of the expansion, iu pi t. 
We note that equation (2) is also nearly identical to the 
expansion potential defined by Samir et al. [1 983 ] and thus 
the two formalisms (charge neutral and nonneutral) provide 
nearly identical downstream potential structures. The vari- 
able Az is the vertical distance from the wake flank into the 
void. At f = 0, the plasma void is formed and the plasma 
discontinuity is in the form of a perfect step function located 
at z — 0 (exactly along the wake flank). As time evolves, the 
plasma expands along z’’, moving Az into the void. For a time 
stationary wake, the time f corresponds to a specific distance 
that a solar wind plasma fluid element has conveeted 
downstream from the obstacle, f - Axw pi /V sw with this 
downstream distance, Ax, being the horizontal distance from 
the obstacle edge (see Figure 1). As such, die ambipolar 
potential inside a crater for nominal solar wind conditions 
(n eo - 5/cm’, T — 10 eV, 1 5 m Debye length, Vsw ~ 400 km/s, 
and w p j — 3000/s) is 

e^/fc T e = ~9.7(Az/A x) - 1 (3) 

For example, consider solar wind flowing horizontally at the 
terminator or polar region over an obstruction of 1 km. in 
height. At distance of 500 m behind the obstacle on the 


leeward side of the obstruction, the value of Az/A x is 2 
making die local ambipolar potential from equation (3) at 
<p A - — 20 kTj/e or ~~200 V relative to the solar wind plasma 
potential. The local electron density' is 

n e /n ra —exp/evw/kT*) (4) 

and is reduced by 1 0‘ J compared to solar wind values. At a 
distance of 1 km behind the obstacle, the local electron 
density is reduced by only t O 5 compared to solar wind values, 
indicating that the solar wind plasma is further filled in at this 
distance. 

[ 17 ] Figure 4 shows the ambipolar potential within 
Shoemaker crater (i.e., our test profile), assuming the solar 
wind flow is horizontal and out of the wrest. The initial solar 
wind conditions that define the properties of the wake 
expansion are n e „ = 5/cm’, 10 e V temperature, and V sw ~ 
400 km/s, implying a 15 m Debye length and u> pi ~ 3000/s, 
This case was run with the solar wind at a 0° elevation angle 
(horizontal flow) relative to the overall surface. The starting 
point for the expansion into the void is the relative maximum 
in the surface located at the x = 18.6 km. As indicated by the 
dark regions, the ambipolar potential is largest along the 
leeward facing wall. Figure 5 shows the ambipolar potential 
values at the floor of the crater (topographic profile also in- 
cluded). Note that at the leeward edge of the crater (Figure 5, 
left), the near-surface ambipolar potentials can drop as low as 
-90 V. This ambipolar potential acts to retard thennal elec- 
tron flow and force ions to catch up to the faster electrons. 


4. Surface Potential 

[is] The surface potential along the crater floor is defined 
by the currents that reach the surface after passing through 
the ambipolar potential. In essence, the ambipolar potential 
“filters” the solar wind reaching the surface. A second near- 
surface potential (the sheath potential) then forms to create 
current balance at the surface. We now address the forma- 
tion of this second potential. 

[is] To define the ambipolar currents that reach any point 
within a PSC, we apply a Samir et al. [1983] expansion of 
Maxwellian plasma with n c = n ( — n co expie r c v kT,.), where 
n C() is the density of the ambient free-flowing solar wind 
away from the obstacle and we initially assume a neutral 
plasma. The ion expansion velocity is (Az/f) + C. where C s 
is the ion sound speed (=A D a,pi) and as we demonstrate will 
have a specific flow angle at the surface, defined as ffu ow . 
Samir et al. [1983] and Crow et al. [1975] define the evolu- 
tion of the wake by a time after the formation of the void, t. 
In an assumed time stationary flow, this time can be re- 
placed by the Ax/V sw , where Ax is the horizontal distance 
from the obstacle edge (see Figure 1) and V sw is the 
downstream ion flow velocity [Farrell et al, 1998]. The 
near-surface plasma currents in the ambipolar region for 
near-horizontal solar wind flow are thus described as 

Je ~ n«o ev tt exp(ep A /kT 4 )(! - 4<r) (5a) 

hr ~ --n M> e(/(Az/Ax)V m . + c.) + V sw Sin(i?sw)) exp(e^ A /kT,) 

(5b) 

Jo, ~ n M e V s ,, Cos (<?**) exp(ey> /kT *) (5c) 
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where <p A *= <p A (Ax, Az) is the ambipolar potential calcu- 
lated at the floor of the crater (shown Figure 5) that has a 
gradient in both x and z. As a consequence, the ambipolar E 
field also is a function of both variables. Variable v to is the 
electron thermal velocity which varies as T<!''~. The angle # sw 
is the initial elevation angle of the solar wind relative to the 
x direction, which is assumed to remain relatively small for 
nominal solar wind near the south pole. We assume the 
electron current incident at the surface releases secondary 
electrons with an efficiency <w between 0.2 and 0.8 
[Halekas et al, 2002, 2009a], The geometry defining Az 
and Ax is illustrated in Figure 1. 

[ 20 ] The thermal electrons are assumed to be quasi- 
isotropic, However, the ions form a well-defined and di- 
rected beam that develops a vertically downward component. 
v z , at the expansion speed of (Az/Ax) V sw + C s . Now 1 con- 
sider a solar wind flow that initially is oriented along the x 
direction (0 SW * 0°). As the plasma propagates downstream 
from the obstacle, plasma is drawn into the void region, with 
a rarefaction wave propagating upward into the solar wind 
and plasma expansion downward into the ambipolar region 
[e.g., Samir et al., 1983, Figure 11]. In this ambipolar 
expansion region, the ions flow inward with a perturbed 
deflection angle, 0 nnvi = arctan(J,/lj X ) ~ arctan(( Az/Ax) + 
(C s /V sw )). Note that for distances dose to the obstructing 
source, Az/Ax > > 1 , the ion flow diverts substantially from 


its initially horizontal direction. For example, the depth in the 
zdirection, Az, in Shoemaker crater is -5 km at the floor. At a 
distance of 5 km from the obstruction edge, the ion flow' di- 
verts into the crater at an angle of ~3‘0° due to the strong 
ambipolar fields. In the middle of the crater at 50 km, the 
deflection angle is closer to 1 0° and for large distances from 
the obstruction edge, the angle reaches an asymptotic value 
near 6° from the horizontal (i.e., the ion sonic Mach cone 
angle). 

[ 21 ] As evident in the equations, the electron are repelled 
by the negative ambipolar potential, and are limited to 
reaching the surface by a repulsive force varying as exp 
(eyJA/Me) (factor in equation (5a)). In contrast, the cold ions 
are accelerated by the wake potential into the crater. The 
surface potential is then calculated assuming current balance 
at the surface in a thin sheath planar geometry: 

•L exp{ vp.,/ k.T c ) JJifSin($n 0W f/topo) -- 9 (fo 

where J e and J, are from equation (5) and # topo is the ele- 
vation angle of the local surface. 

[22] Table 1 lists the current densities and surface poten- 
tial as a function of variable x - the location along the crater 
floor in Figure 5. Note that j, can drop below Kf 10 A/m 2 at 
the bottom of the crater (near the leeward edge). As also 


Table 1. Electrical Environment Characteristics as a Function of Location Along Crater Floor'’ 





x (km) 



18.6 

20 

30 

40 

50 

Ji (A/m 2 ) 

4.6 * 10" 51 

7.7 x 10 i! 

I..0 x iO' 9 

5.4 + 10 10 

6.7 * 10 2 

S' A (V) 

-74 

--63 

-41 

-29 

-23 

•A (V) 

-46 

-22 

-18 

-37 

-18 

ipo” (Vj 

+6; +45 

+6; +45 

+6; +45 

+6; +45 

+6; +45 

l/e Dissipation time {$} 

60 

40 

2.8 

5.6 

0.4 


“The variable x is the horizontal distance across Shoemaker crater shown in Figure 3. The initial solar wind conditions that define the properties of the 
wake expansion, are n,„ «■ 5/errr , 10 eV temperature, Vsw ~ 400 key's, and these make a ! 5 m Debye length arid .ic,,, - 3000 Is, and a 0“ elevation, angle of 
the solar wind coming out of the west, identical conditions to those ran for Figures 4, 6, and 7a, The crater edge is x - ! 8.6 km location. 

’’Here -6 V for object with secondary emission Air'~ 9.5 and 445 y for object with secondary emission — 0.99. 
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Figure 6. The total potential (surface + ambipolar) relative to the solar wind along the crater floor for 
conditions applied to Figure 4 and Table I. The values near -180 V are associated with crater floor 
locations where ion flow is inaccessible; the locations are within the electron cloud. We assume upward 
dust currents thus provide closure for current balance. 


indicated in Table I , the overall potential is the largest at the 
leeward edge of the crater: the two electron retarding po- 
tentials (ip A + r'J gives the total potential at the surface 
relative to the solar wind close to ~ - 90 V near x = 18 km. 
Thus, a solar wind electron has to be energetic enough to 
pass through both negative potentials <p A and <p s to reach 
the surface - which greatly limits the local environmental 
currents at the floor of the PSC. The density of electrons 
is thus n e ~ n eo exp(e(v5 A + ip s )/kT e ) ~ 5 x 10 -4 /cm 3 : a 
nearly 10,000 to 1 electron rejection from the original 
solar wind electron flow overtop the crater assuming a 
Maxwellian plasma distribution. 

[ 23 ] While the electrons thermally advance into the crater 
at the thermal speed, the ions have a very specific and di- 
reeled flow. As a consequence, there are distinct locations 
along the crater floor where ion currents simply do not flow 
onto the surface and the regions are thus immersed in the 
Crow et al. [1975] electron cloud. Figure 6 indicates that 
locations like the leeward crater wall and local maximum 
within the crater itself (like near the x = 45 km location) are 
possible topography features where ion flow may become 
compromised. Unfortunately, current balance at these sur- 
face locations is never reached: dQ/dt = J e exp(eip s /'kT c ) < 0 
and the simple solutions are a continuous accumulation of 
negative charge until the potential goes to negative infinity. 
This is physically unrealistic and thus we infer the presence 
of some other “current of last resort” that gets established to 
provide closure for current balance. 

[24] An obvious (but nonunique) choice of currents is that 
from dust grains repulsed from the strong negatively 
charged surface immersed in the electron cloud. The LEAM 
impact detector operating at the Apollo 17 Taurus-Littrow 
valley landing site measured greatly increased dust activity 
at lunar terminator crossings which have been interpreted as 
dust impacts at tens of m/s that were suggested to be accel- 
erated by terminator E fields [Berg et al, 1976; Berg and 


Perkins, 1979], The flow' of these grains is greatest toward 
the night side directions, implying that the solar wind leeward 
edges of the surrounding mountains became strongly nega- 
tively charged (see Figure 2). LEAM detected an impact 
about every two minutes and to be detected the grain had to 
contain about 1(T 12 C of charge [Berg et al . , 1976], LEAM 
was not initially designed to detect these slower moving 
grains and there is some ambiguity' in the determination due 
to the dissipation time of the LEAM pulse height capacitor 
sensor [Berg and Perkins, 1979]. As such, the detection flux 
appears to be correct but the detailed physical properties of 
the grains (speed and charge) remains of some debate. 

[ 25 ] If we assume the quoted Berg et al. [1976] results, 
then the current density for these grains is j Do ~ A" 1 Q d /r ~ 
1 (T 14 A/m 2 , where A is assumed to be the size of the detecting 
system (0.01 nr), Q d is 10 -12 C and r is the dust incidence 
time of ~100 s. We now incorporate this dust current in 
regions where ions are not directly incident. We might also 
suggest that there is some dependency of released dust on 
surface potential (J D = J Do exp (#/>)) but to assume some 
value of scaling potentia l, < 2 %, would make the situation even 
more contrived. Hence, we leave a more advanced synopsis 
of the dust reaction for future work. 

[ 26 ] Figure 6 shows the total (surface and ambipolar) 
potentials relative to the solar wind along the topographic 
profile of Shoemaker crater. We assume the solar wind 
passing over the top of the crater has an initial n et , = 5/cm \ 
T = 10 eV, V sw ~ 400 tan/s, 7 ^ ~ 0.3 and westward flow 
with an. elevation angle, $ sw - 0° (initial ail horizontal flow). 
Note that the surface becomes strongly negative along the 
leeward edge of the crater, dropping to values near -100 V. 
Each of these solutions (values >-120 V) has ions present at 
the surface location, thus allowing current balance and a 
unique solution for the surface potential. The local topo- 
graphic gradient angle and ion flow angle are included in 
the solution, and account for the large variation in surface 
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Figure 7. A 54 km * 54 kin regional view of the surface potential at Shackleton crater (bottom left) and 
the southern portion of Shoemaker crater (top right). The westward flowing solar wind elevation angle 
and solar illumination angle are 0°, 2°, and 10°, respectively. The 10° case is provided as an illustration 
of extreme angles, but solar illumination at such angles never really occurs near the south lunar polar 
regions presented. White regions represent total potentials at the surface of >0 V and dark regions are 
<-50 V. Note that strong negative potentials exist within the shadowed craters, especially along the 
westem/leeward faces. 


potentials along the crater floor. At the crater floor, local 
surfaces that have normals oriented into the ion flow will 
have a decreasing negative potential (larger Jj), while surfaces 
with normals directed away from the ion flow will have larger 
negative potentials (smaller Jj). 

[ 2 ?] However, there are clearly locations where the ions 
do not flow directly to the surface. One such location is near 
the bottom of the leeward -Facing crater wall (x. = 18 to 
20 km in Figure 6). In these cases, we assume upward neg- 
ative dust currents replace downward ion currents, making 
the total potential drop near -180 V (the crosses located near 
the bottom of Figure 6). In this electron cloud region, the 
total potential at the surface relative to the solar wind is 
e<p s - ey? A ( z) - kT e In (J e (1 - £««)/Jdo) and is dependent on 
the assumed dust current (or some other breakdown current). 
Since J e ~ n^,, e v le exp(e^> A /kT e ), one then finds that 
&p s ~ -kT e ln(n co e v te (l - StstfJoo) which is quasi constant. 
Obviously, this Jdo current of last resort applied herein is not 
unique, and we leave open the possibility of some other 
breakdown current that has yet to be identified. 

5. Regional Views 

[ 28 ] The total potential at the surface can be calculated 
across an entire region, taking the steps above and repeating 
them for each topographic structure. Figure 7 shows a 
54 km x 54 km region containing Shackleton crater and the 
southern rim of Shoemaker crater (see boxed area indicated 
in Figure 3). We note that radar signals did not extend into 
the floor of the smaller Shackleton crater and all elevations 
were set to a constant value of 0 km relative to the mean lunar 
radius. This floor value is still well below the crater rim and 
gives rise to ambipolar expansion effects. In this process, 
relative maxima were identified. On the forward, sunlit side 
of the structure, both solar wind and photoemission currents 
dominate and the surface potential is easily derived from, 
current balance. For the sake of demonstration we have 
aligned the solar illumination and solar wind directions. 
However, the code can consider separate radiation and par- 


ticle elevation angles. On the leeward side of the structure, 
the location of the surface relative to the obstruction maxi- 
mum height is determined, and both ambipolar and surface 
potentials are derived using equations (3), (5), and (6). 

[ 29 ] Figure 7 shows the total surface potential relative to 
the solar wind (<p A + J under a solar illumination and solar 
wind elevation angle # sw of 0° (Figure 7, left), 2° (Figure 7, 
middle), and 10° (Figure 7, right) above the western hori- 
zon. We note that the case for 10° is in fact fictitious and 
that solar illumination never extends to that elevation above 
the horizon at this polar location. However, the case is used 
to demonstrate the effect of sun rise on the region and would 
be applicable for a similar crater at midlatitudes crossing the 
terminator. The initial solar wind conditions that define the 
properties of the wake expansion are n ra - 5/cm 3 , 10 eV 
temperature, S c ff — 0.3, Vsw ~ 400 km/s, and these corre- 
spond to a 15 m Debye length and nj pj ~ 3000/s for use in 
equations (3>~(6). The dark regions correspond to potentials 
at values below -50 V. Note that negative potentials are 
located on the leeward sides of orographic structures and 
inside the leeward/westem faces of craters, At 10° elevation 
angle, the region has a larger portion of the surface exposed 
to direct solar wind and solar radiative flux, and hence the 
area of negative potential is correspondingly reduced. 

[so] Figure 8 shows the same region calculated with 0 = 2° 
but for a warmer plasma consistent with that from a solar 
coronal mass ejection (CME). Specifically, a CME consists of 
an explosive release of driver gas from the sun, As the CME 
driver gas propagates toward the Earth at supersonic speeds, a 
bow shock and plasma sheath (similar in nature to a magne- 
tosheath) develops forward of the CME driver gas, The 
ambient low energy plasma in the sheath region consists of 
densities nearly 4 times the ambient solar wind, warmer 
plasma at many 10 s of eV, and faster flowing plasma at 
600-800 km/s. For figure 8, we run the case for ambient 
densities near 20 et/em 3 , plasma temperatures of 40 eV, 
and flow speeds at 600 km/s. Since secondary,' emissions 
vary roughly as T,, in this regime, we also assume a higher 
surface secondary electron coefficient of 6 si r ~ 0.8. Note 
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Figure 8. The same region as Figure ? but now with solar 
wind conditions during a passing CME (n^ = 20 et/cm 3 , 
T ~ 40 eV, and Vsw ~ 600 km/s, solar wind elevation angle 
of 2°, and 6 ea ~ 0.8). Black regions now indicate surface 
potentials <-150 V. 


that as the solar wind becomes warmer, the surface poten- 
tials behind obstructed regions become strongly negative. 
Figure 8 shows that most of the region has a surface potential 
below -125 V. 

[ri] Figure 9 shows a profile along the floor of Shoe- 
maker crater during a passage of a CME. Figure 9 is anal- 
ogous to Figure 6, only now the westward flowing solar 
wind is moving at 600 km/’s, with a temperature of 40 eV 
and density of 20 el/cm 3 . Note that the region of largest 
potential drop is located along the leeward flank of the 
crater, with surface potentials about 3 times greater during 
the CME sheath passage as compared to nominal solar wind 
conditions shown in Figure 6. 

[ 32 ] We note that Figures 8 and 9 are associated with 
changes in the “cold” plasma during a solar storm/CME 


direct passage to the Moon. We do not consider herein the 
effect from solar energetic particles (SEP) events that are 
associated with a solar flare/CME formation. Such energetic 
particles propagating at relativistic speeds arrive at the 
Moon days before the CME driver gas and sheath modeled 
herein. These energetic particles also have the capability to 
create anomalous surface charging to large negative values, 
as has been reported by Halekas et a!. [2007, 2009b]. 
Halekas et al.'s [2009b] recent modeling of the measured 
particle spectrum from low to high energies should enable 
the additional contribution of such high-energy currents 
from SEPs to be included in our model and this will be 
presented in future efforts. We only consider the cold 
plasma variations at this time, to focus on plasma expansion 
in crater-formed miniwakes, 

[ 33 ] As another example, we include a regional view of 
the highlands located near 84°S, 40 C E (Figure 10). These 
mountains extend 6 km in height but are modified by the 
impact basin associate with Nobile crater - giving rise to 
stark contrasts in topography. For the model, we assume 
the solar wind elevation angle is #sw = 1° above of the 
local horizon, flowing directly out of the east. The dark 
region represent locations where the surface potential is 
<55 V and light regions represent locations where the 
surface potential is >0 V. Note that the highlands do in- 
deed act as orographic structures that create ambipolar 
fields and large, negative surface potentials within the 
crater regions. 

6. Objects on the Crater Floor 

[ 34 ] While the surface charges in an attempt to balance the 
upward and downward current flow's, any object placed on 
the surface will become charged in these same environ- 
mental currents. However, because the object has a differing 
secondary electron character and sheath character {large 
compared to object size) the potential of the object, 

will differ from that of the surface. Specifically, the sur- 
face will charge to obtain current balance, and these local 
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Figure 9, The same surface topographic region as Figure 6, now with similar storm conditions as Figure 8. 
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f igure 10. The potential in the 60 km x 60 km region at 
the edge of the highlands located at 84°S, 40°E. This region 
includes an abrupt topographic change defined by the north- 
east rim of Nobile crater. In this case, the solar wind is flow- 
ing out of the east at a 1° elevation angle, with n« D - 5/cm 3 , 
10 eV temperature, Vsw ~ 400 km/s. Note that the dark 
regions are indicative of surface potentials <-55 V while 
white regions are potentials >0 V. The highland surfaces 
facing the sun/solar wind obtain positive potentials while 
the Nobile crater floor develops a strong negative potential . 

environmental currents in the sheath that form along the 
crater floor are 

~ t»co e v Se expi'e^r/kTcK I - 4nj ~ Ji-eov (?) 

where the total surface/solar wind potential drop is <px — ip s + 
ip A . To complete the calculation, we assume that the object is 
not located in the electron cloud and has direct access to the 
local ion beam. We also assume that the object can redis- 
tribute its charge quickly over its surface. Given these 
assumptions, the object is immersed in these environmental 
currents and will become charged as it attempts to balance the 
electron and ion flux incident at its surface. This static current 
balance at the object surface can be expressed as 

JUnv (l - C?) exp(e^„bj/kT e ) - V e „ v = 0 (8) 

The ion current from the incident wake beams is defined as 
J i ..... env ~ e^obj/mv 2 ] where v is the beam velocity. However, 
for a beam energy near 1 keV. <np a b j/tnv~ « 1 and thus has 
a negligible contribution. The object’s secondary electron 
emission is included in the factor $$$. We assume the object 
is small compared to the local Debye sheath and thus the 
object is acted upon by the local currents at the crater floor. 
As indicated in equation (7), electron and ion currents are 
nearly equal at the surface as forced by the surface potential, 
if the object is larger than the sheath, then the environmental 
currents from the ambipolar region are directly incident (i.e., 
those described in equation (5)). These ambipolar currents 
are no longer balanced. We do not consider this case, but 
choose to leave the nomenclature for J e . cm and J;. e „ v in 


equation (8) for completeness even though they approxi- 
mately cancel for this application. Table 1 indicates the 
object surface potential difference assuming two differing 
secondary' electron emission characteristics: 6 ° - 0.99 and 
( 'cff “ 0-5. We note that the more emissive object charges to a 
larger potential than the lunar surface. For poor secondary 
electron emitters, like many insulators, some concern has to 
be given to the redistribution of charge over the object. If that 
is slow, then the object may differentially charge based on 
access to tire ion flow. A lull 3D model of an object is not 
considered herein, but certainly worthy of consideration for 
future efforts. 

[.is] As an object (astronaut, rover) moves over the reg- 
olith, it is anticipated to collect charge via frictional/contact 
electrification effects between the object and regolith. This 
dynamic phenomenon is commonly called triboeharging. 
The charge exchange is associated with the differing contact 
potentials (which are related to the work function) between 
materials [Desch and Cazzi, 2000], Typically, if a metallic 
material comes in contact with an insulating material, the 
metal: will gain a net charge negative while the insulator will 
obtain a net charge positive. The charge amount for smaller 
objects is a direct function of the triboelectric potential 
difference and related to size differences (relative capaci- 
tance) of the object. 

[36] It was recently shown that an object on the lunar 
surface can be modeled via an equivalent circuit consisting 
primarily of a capacitor connected to both the free charge in 
the regolith and the plasma [Farrell et al., 2008b’J. As an 
object roves over the surface, it will become tribocharged. 
However, since the lunar regolith is a poor conductor, the 
object is most easily grounded via the conductive plasma. 
The electrical dissipation time to remove excess tribocharge 
isT~(C kT/eJenvA), where G is the object capacitance, J rav is 
the environmental current of opposite polarity to the object 
charge polarity' (positively charged object draws in environ- 
mental electrons, negatively charged object draws in envi- 
ronmental ions), T is the temperature of the species drawn to 
the object, and A is the current-collecting area of the object. 

[37] Within the very cold shadowed regions, the regolith 
conductivity is very low, ~1CT 14 S/m [Carrier et al., 1991] 
and electrical dissipation into the regolith is very slow. 
However, as we demonstrate in Table 1, local plasma cur- 
rents also become choked off behind obstacles due to the 
passing flow and both the ambipolar and surface potentials. 
As such, electrical dissipation times increase greatly. 

[3«] Because triboeharging is dynamic, the electrical 
system does not. achieve current balance when triboebar- 
ging. Hence, a charging equation is solved: 

G erf-Cj ,■ dt L-ckv A ( i /.,/ j t-Xp { ty.'obj / k i e ) •- J|-. e r,v A 

C -peb) / 03 (9; 

where vy*) is set equal to an initial value ip c from tri- 
bochargitsg. The »? obl /r R term is included in the formalism 
and represents the dissipation of charge into the lunar reg- 
olith. r R is the dissipation time defined by the regolith 
permittivity and conductivity as ~ 1000 s within cold 
PSC [Carrier et al., 1991], 

[ 39 ] Figure 1 1 shows the electrical dissipation at a loca- 
tion of x = 1 8.6 km, 20 kin, and 40 km along the crater floor 
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Figure 11. The electrical dissipation of an object of I rn 2 surface area placed within Shoemaker crater at 
a location of x = 1 8.6 km (near westward wall), 20 km, and 40 km in Figure 6. The solar wind parameters 
are those applied to Figure 6. Table 1 provides further details on dissipation times and equilibrium values. 
The object is initially tribocharged to -50 V relative to the surrounding surface and has a secondary 
emission efficiency <5“|l - 0.5. Note that near the leeward flowing side of the crater (at x - 18.6 km), 
electrical dissipation of the -50 V object takes many tens of seconds due to the low levels of ambient 
plasma currents at levels of a few JO - ’ 1 A/m 2 . Thus, environmental currents cannot easily remediate the 
excess charge buildup at such locations. Closer to the middle of the crater (x = 40 km), plasma currents 
are larger and thus dissipate the charge buildup faster (on time scales of seconds). The equilibrium 
potential is that relative to the total surface potential and is typically +6 to +45 V positive relative to 
the surface potential depending upon the specific value of (see Table 1). 


for an object of current collecting area. A = 1 m 2 . Table 1 
also lists these electrical dissipation 1/e times. These values 
are calculated for nominal conditions like that for Figure 7a 
(n c0 = 5/cm 3 , 10 eV temperature, Vsw ~ 400 km/s, a 15 m 
Debye length, w pi = 3000/s, and few =■ 0° directed out of 
the west). Note that at x - 1 8.6 km, for an initial tribocharge 
of -50 V, it takes over 100 s to return to ambient potential 
values defined by current balance (dcs/dt = 0). The local 
currents at 40 pA/nr are so low that there is little envi- 
ronmental charge to be drawn at the western edge of the 
crater. In contrast, further eastward of this wall, plasma flow 
is not as restricted, allowing more environmental current to 
be diverted into the crater floor via the ambipolar potential. 
As a consequence, dissipation times progressively decrease 
with distance from the leeward-facing wall (Table 1). 
Hence, at x = 40 km, the local currents are ~0.5 tiA/m 2 
and the corresponding dissipation time is on the order of 
5 s. Note that these dissipation times are far slower than 
those in directly sunlit regions where local currents are at 
micro-Amp/m 2 levels and dissipation occurs on time scales 
of milliseconds [Farrell et ai, 2008bj. Also note that these 
times vary inversely with effective radius: an object 10 times 
smaller will take 10 times longer to dissipate charge. Thus, 
a rover wheel may take hundreds of seconds to reach 
equilibrium levels. 

[40] The dissipation times calculated in Table 1 are for 
locations where ions are fully capable of flowing. However, 
the situation may become very complicated in locations 


without ions in the electron cloud where the ambipolar 
expansion processes are charge nonneutral [ Crow et al, 
1975], Consider an object roving in a region where there is 
an electron doud/void of ions. Further, assume this object 
obtains progressively larger negative tribocharge from the 
very insulating regolith. As the object moves over the surface 
and charges negative, the overlying regolith just roved over 
will tribocharge to a more positive potential than that of the 
object. The strongly negatively charged object will repel 
local electron currents, but will attract objects of differing 
potential: like that very dust just roved over. The local dust of 
lower potential thus becomes the “current of last resort" to 
remediate the object’s negative charge buildup. This effect 
may lead to increased dust clinging in electron cloud, regions. 
This regolith attachment is consistent with the picture that 
small charged particulates become a significant current 
source to stem hypercharging by the object. We describe the 
process of dissipat ion in an elec tron cloud herein (for the first 
time-), recognizing that modeling of such an effect is for a 
future effort. 


7. Plasma Expansion and Its Affect on 
Accumulated Volatiles Within PSCs 

[4 1 ] Without the inclusion of ambipolar plasma wake 
expansion processes, solar wind ions would simply pass 
over the polar crater and not easily move into the deep void 
regions. As indicated in Figure 1, ambipolar expansion E 


H of 14 




E03004 


FARRELL ET AL.: MIMWAKES IN LUNAR CRATERS 


E03004 


Table 2. Ion Sputtering Loss Rate of Ice Within Shoemaker Crater for the Electrical Environment Shown in Figure 6 and Table 1* 





x (km) 




18.6 

20 

30 

40 

51 

) 

.1, (A/nT) 

4.6 * 10 i! 

7.7 x 10 !l 

1.0 * 10 " 

5.4 x 10 10 

6.7 x 

10 

", IV. 

1.1 x 10 4 

2.0 x 10 4 

3.0 x 10 3 

1.6 ^ 10 3 

2.0 x 

10”' 

Ion Energy (keV) 

1.27 

1.14 

0.95 

0.88 

0.85 


Y, Yield for Ice Sputtering^ (molecule/incident ion) 

(1.8 

0.8 

0.8 

0.8 

0.8 


F. Ion Flux (I/nri s) 

2.9 x [O' 

4.8 x 10 8 

6.5 x K) 9 

3.3 X 10" 

4.2 x 

1G 16 

<$> - YF 5 Sputtering Ratefl/nT s) 

2.3 x )0 8 

3.8 « 10* 

5.2 x |(f 

2.6 x 10 9 

3.3 X 

I0 !0 

Loss time for 1 f.im layer of ice (kyi) 

4628 

2801 

205 

409 

32 



'Here tij/n w is the ratio of wake ions to the free-flowing solar wind density of 5/cm 5 . The crater edge is x ~ 18.6 km location. Yields are for normal 
incidence. 

"Johnson [1990, Figure 3.22b). 


fields deflect the solar wind inward into the initial void re- 
gions, allowing ions to be incident or strike at the crater 
floor (see Jj, Table I). Such ions may be sources of sput- 
tering which may act to erode local volatiles (e.g,, water ice) 
that have accumulated at the surface of the PSC. 

[42] Table 2 shows the effects of ion sputtering from an 
icy surface on the crater floor, assuming the deflection of 
solar wind protons as modeled in Figure 6 and Table 1 
(n co = 5/cm 3 , 1 0 eV temperature, V svv ~ 400 km/s, 5 dr ~ 0.3, 
6W = 0°). The ice sputtering yields for H + of ~0.8 molecule 
per incident ion are from Figure 3.22b of Johnson [1900], 
Note that sputtering loss of water is greatest in the central 
crater region (x = 50 km) where inflowing ion densities are 
relatively high (at -2% solar wind levels). 

[43] Table 2 also show's the loss rate of a polar ice layer in 
units of kiloyears per micron of water layer thickness. Note 
that in the central region, a water layer is eroding — 120 times 
faster than it is at the leeward edge of the crater. However, 
we note that the leeward edge of the crater changes position 
over the course of a lunation: thus the ion influx varies from 
small to large as the wall surface norma! rotates from anti- 
sunward to sunward directed. Solar illumination will also 
vary at the rim of the crater as a function of lunation as well, 
possibly exposing the rim region to direct sunlight and more 
positive surface potentials (e.g., see the x = 70 km region in 
Figure 6). Thus, the low erosion rates along the crater wall 
are only realized for a fraction of each lunation. 

[44] While we model the relatively large Shoemaker cra- 
ter, we note that small but deep polar craters may experience 
the least overall ion erosion, as suggested by Table 2. For 
example, consider a case of a 5 km deep crater that extends 
only 20 km in diameter (x = 40 km position in Table 2). In 
this ease, the water loss rate in the crater never drops below 
0.4 My per micron anywhere along the crater floor, and 
erosion via ion sputtering will take a relatively long time. In 
larger craters greater than 30 km (x = 50 km position in 
Table!) overall ion inflow to the floor is greater resulting 
in substantial ion/surface incidence which speeds up erosion. 

[45] However, the situation may be different if the sub- 
strate layer is the regolith itself (e.g., Si0 2 with no 0 %'erlying 
ice layer). In this case, the sputtering yields are low, on the 
order of -0.03 molecules per incident ion for ~1 keV H r 
incidence [see Johnson, 1990, equation (3.15), Figure 3.16b, 
and Figure 3.17]. These values are about a factor of 50- 
100 times lower than the yield of ice. Such values suggest 
that any solar wind ion incident on the polar regolith will 
not release material but instead become mass accumulated. 
As such, with a nonice lunar regolith surface, ion deflec- 


tion into the crater becomes a source of l-V at the surface 
[Feldman ex al., 2001; Crider and Vondrak , 2003], 

8 . Conclusion 

[ 46 ] The case is made that large-scale structures along the 
lunar terminator and polar regions can obstruct the solar 
wind (referred to herein as “solar wind orographic effects”), 
creating a wake region in the downstream flows. These 
plasma miniwakes are subject to well-known plasma am- 
bipolar expansion processes that define the electric potential 
structure at the surface of leeward mountain faces 'and within 
polar shadowed craters. On the leeward sides, plasma cur- 
rents are vastly reduced, and ambipolar potentials develop as 
electrons migrate into the shadowed regions ahead of 
slower, massive ions. We apply analytical models [ Crow et 
at , 1975; Samir et at, 1983] (that are also consistent with 
simulations; Farrell et at [2008a]) of the ambipolar process 
to define the near-surface ambipolar potentials and currents 
reaching the leeward side of lunar obstacles in a near- 
horizontal solar wind flow at the lunar polar region. These 
ambipolar-created currents are then used in a second calcu- 
lation to derive the potential directly at the surface where all 
currents are balanced to equal values. We note that there are 
regions along the crater floor where the expansion loses 
neutrality (ion flow cannot reach the surface). At these lo- 
cations, an electron cloud is fanned. The formation of the 
electron cloud is also consistent with kinetic simulations of 
the expansion [see Farrell et at, 2008a, and references 
therein], 

[47] Within lunar craters, the surface potential can be 
derived in most locations. However, those locations lacking 
ion flow (i.e., electron cloud regions) cannot be properly 
balanced without the inclusion of some new downward 
positive/upward negative current. There are a number of 
possible choices for this remediative current, and one ob- 
vious (albeit nonunique) source is the current produced by 
lofted negatively charged dust from regions within the 
electron cloud. Observations from Apollo 17*s surface 
LEAM experiment indicate that dust activity increases at 
terminator crossings where local orographic effects from 
nearby mountains are the largest, Hence, there is Apollo 
observational evidence that we might expect surface-ejected 
dust to provide the needed currents. Using LEAM results (as 
a starting point) we suggest that dust provides an upward 
negative current of >10‘‘ 4 A/rrri. Using this value, we can 
obtain the needed current balance - and find that the surface 
potentials, tp u can become more sharply negative within the 
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electron cloud regions, relative to the nominal solar wind 
flowing over the top of the shadowed region (see Figure 6). 

[4s] We suggest that lifted dust is not simply a “science 
curiosity” hut may provide a critical role in the closure of 
current balance in ambipolar-created electron cloud regions. 
Negatively charged dust is lofted from the hypercharged 
surface covered by the electron cloud in an attempt to re- 
mediate the excessive charge buildup at the surface/doud 
interface. The charged dust becomes the current of last re- 
sort that is expelled in an attempt to achieve current balance 
(attempts to reduce dQ/dt). As such, this picture of ejected/ 
lofted dust fits naturally in the overall picture of the elec- 
trical environment within polar shadowed craters and behind 
obstructing mountains. We would thus predict that LEAM’s 
“dusty sleet” should be anticipated in the polar environment 
that is essentially in a perpetual terminator. Clearly, some 
sort of precursor validating observation is required to con- 
firm this prediction, but there is certainly an analog with the 
equatorial Taurus-Littrow valley where energetic dust at 
terminator crossings was reported. We note that there may 
be other ways to obtain current balance (increased second- 
ary electron emission, field emissions, etc.) and these need 
to be investigated in the future as well. 

[49] Given the overall electrical environment, we can then 
consider the electrical fate of objects placed within the 
PSCs. We find that stationary' objects tend to charge positive 
relative to the surface (still negative relative to the over- 
flowing solar wind), with the largest differential charging 
occurring with objects that have a secondary electron co- 
efficient that has the greatest difference from the lunar 
regolith. If the object has secondary electron emission 
similar to the lunar regolith, the differential potential be- 
tween the object and lunar surface is reduced almost to zero. 

[so] Of concern is electrical dissipation of any tribocharge 
buildup from a moving object, since both lunar surface 
currents and plasma currents are reduced in PSCs. It was 
found that dissipation time scales can become very large 
near the leeward face of PSCs with equilibrium times 
achieved in 1 00 s of seconds in regions of low ion currents 
(leeward crater walls). Within the electron cloud, we again 
demonstrate that plasma currents will likely have difficulty 
obtaining balance at the object. We suggest that the recently 
tribocharged regolith that has just been churned by the ob- 
ject then act as the remediation current of “last resort.” 
However, this current may still not dissipate charge as fast 
as bleeding the object charge directly into the low conduc- 
tivity ground (i,e., directly into the insulating regolith). As 
such, we expect very slow dissipation times ( — 1 000 s) 
within electron cloud regions of PSCs. Certainly, if any 
human object is roving in the leeward edges of the craters, 
special care must be made for avoidance of electrostatic 
discharge hazards. 

[ 51 ] The modeling, based on previous analytical [ Crow 
etal., 1975; Samir el al., 1983; Martha, 1973; Stubbs et 
al., 2006 ] and simulation [Farrell et al. , 1998; Birch and 
Chapman, 2001a, 2001b; Farrell et al, 2008a] studies, ap- 
plies ambipolar and wake expansion in the vertical plane. As 
an example, for a mountain we consider the expansion as a 
function of vertical dimension. However, we do not consider 
expansion about the edges of a feature (like a mountain), in a 
horizontal plane parallel to the surface. A 3-D PIC code could 
model the multidimensional system and should be consider 


for future work. Such a model may also include plasma 
transport from implanted ions on the windward side and their 
transport over the topographic feature into the expansion 
region. This “ion leakage” might also be a current to re- 
mediate large surface potential regions. We do not mode! 
such effects herein. We also assume a Maxwellian plasma 
distribution describes the ambient solar wind and expansion 
region and future efforts could incorporate a kappa distri- 
bution [Halekas et al. , 2005], Ideally, a multidimensional 
kinetic simulation can be used to consider expansion about 
a complex object. However, the results here provide a first- 
order estimate of surface potentials. 

[ 52 ] We also assume that the plasma is flowing uncon- 
strained by local magnetic fields: that there are no surface 
magnetic fields in the PSC region. Such fields (as yet to be 
identified) will alter plasma geometries and currents, with 
the result in most locations being to further reduce the 
currents reaching the surface [Reiff, 1976). As an example, 
as solar wind electrons move close to the surface and into 
increasingly converging magnetic field geometries, elec- 
trons of large pitch angle are reflected at the mirror points 
and thus never reach the surface. As such, there would be an 
additional, electron and ion retardation factor to equations (5) 
and (6), further increasing dissipation times via plasma 
return currents. 

[53] We disconnect the solution of the ambipolar poten- 
tial, i p A , to that of the surface potential, and this is 
allowable due to the development of the Debye sheath at the 
surface. The sheath plasma is specifically formed to cancel 
the electrical forces from the surface and as such the surface 
charge does not affect the formation of the ambipolar 
expansion (except at regions where Az< A|>). This expansion 
proceeds with minimal influence from the surface below. As 
such, the two processes can be calculated independently and 
solutions merged piecemeal. This situation is unlike an infi- 
nitely extended surface in free space, where any surface 
charge alters potentials out to infinity. In the lunar case, the 
surface is shielded by a conductive (plasma) medium that 
polarizes charge near the surface (the sheath charge) to 
exactly cancel the surface charge beyond a few Debye lengths. 
The potential from this surface goes to zero at infinity. 

[ 54 ] Secondary electrons have been lumped into a bulk 
parameter, S e a, but are of special interest. Specifically, they 
are released from the surface in low concentrations: the 
densities are low enough (n scc ~ n e „ exp(e^A/kT e ) for 

at the surface) that they should not significantly modify die 
ambipolar potential structure at altitudes above the surface. 
We assume that this is the case, but note that even weak 
beams can affect local potentials in nonlinear ways and such 
should be examined in more detail. The beams are accel- 
erated through the sheath and ambipolar potentials to form 
well-defined beams at very high altitudes, like those de- 
tected by Lunar Prospector [Halekas et al. , 2002, 2009a]. As 
such, these emissions become an on-orbit remote sensing 
tool to determine the underlying near-surface potentials. 

[ 55 ] An assertion of this work is that dust is most likely 
and easily lofted from crater floors in regions where the 
surface potential is the largest. Hence, a simple way to 
sample resources within the lunar craters is to perform dust 
collection from orbit — and let the crater surface electrical 
effects move material from the surface to orbiting altitude. 
In essence, we let the lunar environment provide the “heavy 
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lifting” of material so that a sampling system does not have 
to perform a difficult landing and return liftoff within the 
cold, electrically complicated crater. 

[56] The situation may also be true at a base located near i 
above the PSC. Because of electrical repulsion effects, 
LE AM -like dust containing internal crater resources (i.e., 
water) may be ejected out of the crater floor and onto the 
surrounding surface. This effect will maximize when the 
solar wind is flowing in the direction over the crater and 
the lunar base is downstream. In this case, the surface 
normal of the leeward edge of the crater points toward the 
lunar base - directing material in the direction of the base 
itself. As such, explorers may not have to go into the cold 
craters to determine if resources are present: they may 
simply “pop out” of the crater to them. The question is if the 
amount being ejected is enough to be of any real value to a 
landed exploration team. A precursor environmental mis- 
sion would answer this question. We thus conclude that the 
environment has significant electrical variability: but this 
same electrical variability may be used to the advantage of 
any exploration team. 
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